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Abstract 

We study equilibrium liquid crystal configurations in three-dimensional 
domains, within the continuum Landau-De Gennes theory. We obtain 
explicit bounds for the equilibrium scalar order parameters in terms 
of the temperature and material-dependent constants. We explicitly 
quantify the temperature regimes where the Landau-De Gennes predic- 
tions match and the temperature regimes where the Landau-Dc Gennes 
predictions don't match the probabilistic second-moment definition of 
the Q-tcnsor order parameter. The regime of agreement may be inter- 
preted as the regime of validity of the Landau-De Gennes theory since 
the Landau-De Gennes theory predicts large values of the equilibrium 
scalar order parameters - larger than unity, in the low-temperature 
regime. We discuss a modified Landau-De Gennes energy functional 
which yields physically realistic values of the equilibrium scalar order 
parameters in all temperature regimes. 
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figurations, Landau-De Gennes theory 
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1 Introduction 

Liquid crystals are an intermediate phase of matter between the commonly 
observed solid and liquid phases [6]. In the simplest liquid crystal phase, 
the nematic phase, the constituent rod-like molecules translate freely as in 
a conventional liquid but whilst flowing, tend to align along certain locally 
preferred directions i.e. they exhibit a certain degree of long-range orien- 
tational ordering. Liquid crystals have attracted a lot of interest in recent 
years because of their unique physical properties and continue to do so be- 
cause of their diverse applications [12] and their analogies to other physical 
systems. 
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The mathematical theory of nematic hquid crystals is very rich; for a re- 
view see |18l [T3] . The key ingredient of any mathematical theory for nematic 
liquid crystals is the definition of an order parameter that distinguishes the 
ordered nematic phase from the disordered isotropic liquid phase. Mean- 
field liquid crystal theories, such as the Maier-Saupe theory, describe the 
liquid crystal configuration in terms of a probability distribution function 
tp on the unit sphere. The order-parameter, known as the Q-tensor order 
parameter, is defined in terms of the second moment of ijj [U [16]. This prob- 
abilistic second-moment definition naturally requires Q to be a symmetric, 
traceless 3x3 matrix and imposes certain constraints on its eigenvalues, 
which represent the degree of ordering. The Landau-De Gennes theory, on 
the other hand, is a continuum theory for nematic liquid crystals and does 
not contain any information about either tp or the intermolecular interac- 
tions [HI [16]. The Q-tensor order parameter, within the Landau-De Gennes 
framework, is a symmetric, traceless 3x3 matrix with no a priori con- 
straints on the eigenvalues. The Landau-De Gennes energy functional is a 
nonlinear integral functional of Q and its spatial derivatives and the equi- 
librium, physically observable configurations correspond to either global or 
local minimizers of this energy subject to the imposed boundary conditions. 

A natural question of interest is - do the equilibrium configurations pre- 
dicted by the Landau-De Gennes theory agree with the probabilistic second- 
moment definition of Q? We systematically address this question in this pa- 
per. We obtain explicit bounds for the scalar order parameters of global en- 
ergy minimizers, referred to as equilibrium scalar order parameters, in terms 
of the temperature and the material-dependent constants. These bounds 
quantify (to some extent) the competing effects of the different terms in the 
Landau-De Gennes energy density. Further, these bounds are compared to 
the probabilistic second- moment definition of Q. This allows us to explicitly 
delineate the regions of agreement and the regions of disagreement and we 
find that the Landau-De Gennes predictions don't match the probabilistic 
second-moment definition in the low-temperature regime. In particular, the 
equilibrium scalar order parameters, within the Landau-De Gennes frame- 
work, can take physically unrealistic values (larger than unity) in the low- 
temperature regime. Our results largely depend on the use of maximum 
principle type of arguments for nonlinear elliptic systems of partial differ- 
ential equations and they can be readily extended to more general energy 
functionals than the ones considered in this paper. 

The derivation of the Landau-De Gennes energy density is valid near 
the isotropic state, close to the nematic-isotropic transition temperature. 
Therefore, it is well-expected that the predictions are physically unrealistic 
in the low temperature regime. However, our results show that the Landau- 
De Gennes predictions fail to be consistent with the probabilistic second- 
moment definition within a sufficiently small neighbourhood of the nematic- 
isotropic transition temperature. In principle, one would want to develop a 
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continuum theory that works for all temperature regimes. In the last part 
of the paper, we briefly outhne a Ginzburg-Landau approach that remedies 
the flawed predictions in the low temperature regimes. We deflne a modifled 
Landau-De Gennes energy functional such that the energy density blows up 
whenever the liquid crystal conflguration violates the constraints imposed 
by the probabilistic second-moment deflnition of Q or equivalently whenever 
the scalar order parameters are physically unrealistic. One deflciency of 
this approach is that it has no apparent connection with the mean-fleld 
microscopic approaches. A different approach has been suggested in [TT] and 
we hope to systematically investigate a microscopic-macroscopic derivation 
of a continuum energy functional in future work [2]. 

The paper is organized as follows. In Section O we review the probabilis- 
tic second- moment deflnition of the order parameter. In Section [3l we study 
equilibrium liquid crystal conflgurations within the continuum Landau-De 
Gennes theory. In Section 3.1, we consider spatially homogeneous cases 
whereas in Sections 3.2 and 3.3, we include spatial inhomogeneities into the 
model and obtain upper bounds for the corresponding equilibrium scalar 
order parameters. These bounds explicitly deflne the domain for the equilib- 
rium scalar order parameters in terms of the temperature and the material- 
dependent parameters. In Section HI we discuss the main results and con- 
clusions of this paper and suggest future research directions. 

2 The Probabilistic Second-Moment Definition 

In this section, we briefly review the probabilistic second-moment deflni- 
tion of the Q-tensor order parameter and the Maier-Saupe mean-fleld liquid 
crystal theory. The interested reader is referred to [6l [TH [T9] for details and 
we present the main points here for completeness. 

Within the simplest microscopic model, the nematic molecules are mod- 
elled by elongated rods where the long molecular axes tend to align along 
certain locally preferred directions [HJ [16]. The state of alignment of the 
nematic molecules is described by a probability distribution function for the 
molecular orientations on the unit sphere, tp : S'^ ^ M"*", since 5^ C is 
the space of all admissible directions. The probability distribution function, 
■0(p), gives the probability of flnding molecules oriented in the direction 
p G S"^ . Then tf; has the following properties - [Si [1] - 



where ([1| accounts for the indistinguishability of the states p and — p on 
the unit sphere. 



V'(P) > peS^ 




(1) 
(2) 
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The macroscopic variables are defined in terms of the moments of ip. 
The first moment vanishes because of the equivalence between antipodal 
points, p = —p. We define the nematic order parameter, the Q-tensor 
order parameter, to be the normalized second moment of the probability 
distribution function as follows (BJ [16]- 

Q = (p ® P - ^l) V'(P) dp (3) 

We refer to ([3]) as the probabilistic second-moment definition of Q in the 
rest of the paper. For an isotropic system, where all directions in space are 
equally likely, the function ^ is a constant i.e. 

V'(P) = ^ VpG52 (4) 

and consequently, Q = 0. On the other hand, for a perfectly aligned system 
where the nematic molecules identically align along a pair of unit-vectors 
(e, — e), the function tp is given by - 

^{p) = ^{Ss2{e,p)+6s2{-e,p)) (5) 

where 6^2 is the Dirac-delta function on S'^ and the corresponding Q-tensor 
is Q = (e(g)e- il). 

It follows directly from ([3]) that Q is a symmetric, traceless 3x3 matrix. 
From the spectral decomposition theorem, we can express Q in terms of a 
triad of orthonormal eigenvectors, {61,62,63}, and corresponding eigenval- 
ues, {Ai, A2, A3}, subject to the tracelessness condition Aj = 0. 

Q = A161 (gi 61 A262 62 -h A3e3 (g) 63 where ^ Aj = 0. (6) 

i 

Nematic liquid crystals are broadly classified into three main families ac- 
cording to the eigenvalue structure of Q. A nematic liquid crystal is called 
isotropic when it has three equal eigenvalues (the tracelessness condition im- 
plies that Q = 0), uniaxial when it has a pair of equal non-zero eigenvalues 
and biaxial when it has three distinct eigenvalues O |16]. The eigenval- 
ues measure the degree of orientational ordering along the corresponding 
eigenvectors and one can verify that the eigenvalues are constrained by the 
following inequalities - 

-\<h= I {v■^^f^{p) dv-\<\, i = 1...3 (7) 
6 J s'z 6 6 

since < (p ■ Gj)^ ^{v) dp < 1. If the eigenvalue Aj = — | (the lower 
bound in ([7])), then the function ip is supported on the great circle perpen- 
dicular to the corresponding eigenvector 6j. On the other hand, if Aj = | 
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(the upper bound in ([7])), then tp is as in ([5]) and the hquid crystal molecules 
line up perfectly along the pair of unit-vectors (ej,— ej). For example, the 
liquid crystal state, (Ai, A2, A3) = (|, — |, — |), is an example of a perfectly 
ordered state along the eigenvector ei and exhibits prolate uniaxial sym- 
metry whereas the liquid crystal state, (Ai,A2,A3) = (— has the 
molecules aligned in the plane orthogonal to ei and exhibits oblate uniaxial 
symmetry [5]. From a physical point of view, the limiting values, Aj = — ^ 
or Aj = |, represent unrealistic configurations. 

The Q-tensor order parameter can be expressed more concisely in terms 
of just a pair of eigenvectors (ei,e2) and a pair of scalar order parameters 
(s,r) as shown below [16j . 

Q = s (^ei ei - ^I^ + r es - ^I^ , (8) 

where s,r are linear combinations of the Aj's given by 

s = Ai - A3 = 2Ai + A2 

r = A2 - A3 = Ai + 2A2. (9) 

The constraints ([7]) directly translate into constraints for the scalar order 
parameters (s,r) in ([9]) and necessarily imply that (s, r) take values inside 
or on the boundary of the physical triangle, T^, illustrated in Figure [H On 
each of the boundary segments of T^, one of the eigenvalues Aj necessarily 
attains the lower bound in ([7]). For example, on the boundary segment 
s + r = 1, we have A3 = — ^. Similarly, every vertex of represents a state 
of perfect alignment along of the eigenvectors of Q. For example, the vertex 
(s, r) = (1,0) represents a state of perfect alignment along the eigenvector 
ei. We call T^, the physical triangle, on the grounds that the scalar order 
parameters are appropriately bounded (less than unity) inside and the 
boundary represents physically unrealistic liquid crystal configurations. 

For definiteness, we can assume a specific ordering of the eigenvalues such 
as A3 < A2 < Ai. Then Ai is necessarily non-negative and A3 is necessarily 
non-positive and the constraints ([7]) require (s, r) to take values inside a 
subset of T^, which is referred to as a fundamental domain Tf defined below 

Tf = {{s, r); < s < 1, < r < min {s, 1 - s}} C T^. (10) 

Analogous remarks apply to the other five possibilities for the ordering of 
the eigenvalues. 

The Maier-Saupe theory is a mean-field theory for uniaxial nematic liquid 
crystals [6l[T7]. The Maier-Saupe free energy has two contributions - 

hism = [ V(p)logV'(p) dp - l-U{T)S^ (11) 
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where U{T) accounts for the intermolecular interactions and is temperature- 
dependent and S is the uniaxial scalar order parameter. A standard mini- 
mization procedure for Ims yields a self-consistent equation for the equilib- 
rium order parameter, S{T), as a function of the temperature. For high tem- 
peratures, the isotropic phase = is the global energy minimizer whereas 
for temperatures below a certain critical temperature Tc, the nematic phase 
is globally stable and the Maier-Saupe theory predicts a first-order nematic- 
isotropic phase transition at the critical temperature Tc- 

r 

I 
I 
I 




(-1,-1) 

I 
I 
I 
I 

Figure 1: The physical triangle T^. The origin (s, r) = 
(0, 0) represents the isotropic liquid state; the dotted lines U = 
{(s, r) G : s = or r = or s = r} \ (0, 0) represent uniaxial states and 
S = \ {[/ U (0, 0)} is the biaxial region. 



3 The Landau— De Gennes Theory 

In this section, we study equilibrium nematic configurations within the con- 
tinuum Landau-De Gennes theory. The Landau-De Gennes theory describes 
the state of a nematic liquid crystal by a macroscopic order parameter - the 
Q-tensor order parameter, which is defined in terms of macroscopic quan- 
tities such as the magnetic susceptibility [H [13]. Within the Landau-De 
Gennes framework, Q is a symmetric, traceless 3x3 matrix with no a priori 
bounds on the eigenvalues; in particular, the eigenvalues are not constrained 
by the inequalities ([7]). 

We work in a three-dimensional setting, take the domain C to be 
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bounded and simply-connected with smooth boundary. Let 5*0 denote the 
space of symmetric, traceless 3x3 matrices 

So = {Qe M3x3; Q^^ = Q^^, Q„„ = 0} . (12) 
The corresponding matrix norm is defined to be [13] 

IQI' = Qa/3Qa/3 a,/3 = 1...3 (13) 

and the Einstein summation convention is used here and elsewhere in the 
paper. We define our admissible space A to be 

A = {Q e W^'^ {n, So) ; Q = Qo on dn} (14) 

where the Sobolev space W^'^ (i7, 5o) is given by [IHl E] 

W^'\n,So) = S^Qe So; y IQI V |VQ|2 < cx)| (15) 

and Qo is a smooth, physically realistic boundary condition in the sense that 
its scalar order parameters, (s, r) in ([U]), are inside the physical triangle T^. 

In the absence of external fields and surface energies, the Landau-De 
Gennes energy functional, Ilg^ is given by 

Ilg [Q] = f fB (Q) + L|VQ|2 dV. (16) 
Jn 

Here is the bulk energy density, L > is a material-dependent elastic 
constant, 

|VQ|^ = Qjj,fcQij,fc i,j,k = l...3 

is the elastic energy density where Qij^k = denote the first partial 

derivatives of Q. 

Comment: We work with the simplest form of the elastic energy density 
- the one-constant elastic energy density in ilb]) . There are more general 
forms of the elastic energy density, see flfij. [^. 

The bulk energy density fs is a scalar function of Q and it dictates the 
preferred liquid crystal phase - isotropic, uniaxial or biaxial. We work with 
the simplest form of fs that allows for a first-order nematic-isotropic phase 
transition. This simplest form of fs is a quartic polynomial in Q as shown 
below 

fB (Q) = ^trQ2 - ^trQ3 + | (trQ^)' with (17) 

trQ^ = Qal3Qa/3, trQ^ = Qa/3Q./3^Q.-/a O, /?, 7 = 1 • • • 3. (18) 

Here b,c > are material-dependent bulk constants, independent of the 
temperature, whereas the parameter a scales linearly with the absolute tem- 
perature and is given by 

a = a{T-T*) (19) 
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where a > and T* is a characteristic hquid crystal temperature [6| I16j. 

The equihbrium, physically observable configurations correspond to global 
or local minimizers of the Landau-De Gennes energy functional, Ilg^ sub- 
ject to the imposed boundary conditions. In what follows, we first consider 
spatially homogeneous cases in Section 13.11 and then study global energy 
minimizers in spatially inhomogeneous cases in Sections 13.21 and 13.31 



3.1 The Bulk Energy Density 

Our first proposition concerns the stationary points of the bulk energy den- 
sity. Proposition [T] is known in the literature [9] and we give an alternative 
proof here for completeness. 

Proposition 1. [Jj The stationary points of the bulk energy density, fs in 



77[), are given by either uniaxial or isotropic Q-tensors of the form 



Q = d (^n (g> n - ^ij (20) 

where d is a scalar order parameter and n is one of the eigenvectors of Q in 
On comparing [2U\) with (0), we see that when n = ei, the parameter 
d = s and the scalar order parameter r = 0. Similarly, when n = e2, the 
parameter d = r and s = whereas when n = 03, the order parameters s,r 
in ^ are equal and are given by s = r = —d. 

Proof. For a symmetric, traceless matrix Q of the form dH), trQ" = Yl\^i A" 
subject to the tracelessness condition so that the bulk energy density fs only 
depends on the eigenvalues Ai, A2 and A3. Then the stationary points of fs 
are given by the stationary points of the function / : ^ M defined by 

/(Ai,A2,A3) = |^A?-^^Af + |f^AA -25^A. (21) 

4=1 i=l \i=l / i=l 

where we have recast fs in terms of the eigenvalues and introduced a La- 
grange multiplier 5 for the tracelessness condition. 

The equilibrium equations are given by a system of three algebraic equa- 
tions 

-Q^aXi- bXj + c I ^ \l \ A, = 2(5, for i = 1 . . . 3, (22) 



9A,; 



\k=l 



together with the tracelessness condition ^ ■ Aj = 0. The system (|22|) is 
equivalent to 



(Ai - A,) 



3 

b[Xi + \j) + cY,>l 

k=l 



1 < i < J < 3. (23) 
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Let {Aj} be a solution of the system (I22p with three distinct eigenvalues 
■^1 / ^2 / A3. We consider equation (p3]) for the pairs (Ai, A2) and (Ai, A3). 
This yields two equations 



MAi + A2) + c^A^ = 

k=l 
3 

6(Ai + A3) + cJ^Ai = (24) 



k=i 



from which we obtain 



6(A2-A3) = 0, (25) 



contradicting our initial hypothesis A2 / A3. We, thus, conclude that a sta- 
tionary point of /b must have at least two equal eigenvalues and therefore 
correspond to either an uniaxial or isotropic liquid crystal state. In partic- 
ular, there are no biaxial stationary points for the particular choice of /b in 
dlTD. □ 



By virtue of Proposition [H it suffices to consider uniaxial Q-tensors of 
the form 

Q = s (g) n - ^I^ ne 

whilst computing the stationary points of /b- For such Q-tensors, /b is a 
quartic polynomial in the uniaxial scalar order parameter s and the station- 
ary points are the roots of the algebraic equation given below 

^ = J- (l8as - 665^ + I2cs^) = 0. (26) 
as 2" 

There are precisely three stationary points; 



b±Vb'^ - 24ac 

s = and s± = — (27) 

4c 



where 



s2 



/b(0) = and /5(s±) = -f (9a-6s±), (28) 

54 

and fB{s-) > /b(s+)- Hence, the global bulk energy minimizer is either the 
isotropic state Q = or the ordered nematic state 

Q = s+ (^e e - (29) 

where e is the eigenvector with the largest eigenvalue. 

A natural question is - for which temperature ranges does the global 
bulk energy minimizer lie inside the physical triangle i.e. for which 



9 



temperature regimes does s+, which is the stable nematic stationary point, 
take values in the physical range 



4c 

One can directly verify that s+ G [0, 1] if and only if 

^(6-2c)<a<^, (30) 

or equivalently, in terms of the absolute temperature T if and only if 
1 62 

— {b-2c) + T*<T<-— + T*. (31) 
6a 24ac 

For the common liquid crystal material MBBA, the values of the character- 
istic bulk constants are given in the literature |15t [T6] 

a = 0.42 X loV/m^ °C, b = 0.64 x lO"^ J/m^, c = 0.35 x 10^ J/m^ 
T* = 45°C Tc = 46°C (32) 

where Tc is the nematic-isotropic transition temperature. We substitute 
these values into §0^ and ^ and find that s+ > 1 for T < 44.52°C i.e. 
s+ moves outside the physical range within a 2°C - neighbourhood of the 
nematic-isotropic transition temperature. 

We recall that there are three characteristic temperatures predicted 
by the quartic form of Jb in ((TTj): (i) a = 0, below which the isotropic 
state loses its stability (ii) the nematic-isotropic transition temperature, 
a = a{Tc - T*) = for which /s(s+) = /s(0) and (iii) a = ^ above 
which the ordered nematic stationary points are no longer defined in ()27p . 
We provide a pictoral representation for the stationary points of /b for 
ease of comparison with T^. We define the bulk triangle, A(T), to be the 
convex hull of the stationary points of /b in the order-parameter (s, r) 

- plane. For -aT* < a < -g, A(r) is an isosceles triangle with its 
vertices at the points {(2|s_|,0), (0,2|s_|), (— 2|s_|, — 2|s_|)} whereas for 

— 1^ < a < A(T) is an isosceles triangle with its vertices at the points 
{(s-(_,0), (0,s+), (— s+,— s+)}. For a > A(T) collapses to the origin 
since s = is the unique critical point in the high-temperature regime. In 
Figures (21 and El we illustrate A(r) for all temperature regimes. 



3.2 The One-constant Elastic Energy Density 

In this section, we study global minimizers of the Landau-De Gennes energy 
functional, Ilg in (fT6]) . 

IlgIQ] = [ /B(Q)+i|VQ|2 dV 
Jn 
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(a) (b) 

Figure 2: (a) The triangle A(r) for + T* < T < T*. The red marked 
points label the stationary points of /b in this temperature regime, (b) The 
triangle A(r) for T < + T*. 




(a) 

Figure 3: The triangle A(T) for T* <T < + T* . The stationary point 
s_ > in this temperature regime, (b) The triangle A(T) shrinks to the 
origin for T > + T* , since the isotropic state (s, r) = (0, 0) is the 
unique stationary point in this temperature regime. 
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and obtain explicit bounds for the equilibrium scalar order parameters, {s, r) 
in ([9|). These bounds quantify the effect of the elastic energy density on the 
bulk energy minima and can also be compared to the probabilistic bounds 
in dZD. 

3.2.1 Existence and regularity of minimizers 

There exists a global minimizer, Q*, of Ilg in the admissible class A, where 
A has been defined in This is a ready consequence of the direct methods 
in the calculus of variations [HIS]. Indeed, one can check that Ilg satisfies 
the following coerciveness estimate [5] 

hclQ] > ac\\Q\\ln,2^n) (33) 
where etc > and the H^^'^-norm, ||Q||vKi.2(rj)5 is given by 



IQIlTyi.^(n) = (^|Q|' + |VQ|2 dV^ 



1/2 



The Landau-De Gennes energy density is convex in the gradient VQ and 
therefore Ilg is weakly lower semicontinuous [7] . The coerciveness and weak 
lower semicontinuity of Ilg guarantee that the infimum energy is actually 
achieved i.e. there exists a Q* G ^ with the property 

IiG[Q*]= inf /lg[Q]. (34) 

The global minimizer Q* is a weak solution of the corresponding Euler- 
Lagrange equations, which is a system of nonlinear elliptic partial differential 
equations as shown below - 

2LAQal3 = aQal3+b QtrQ^^a/? - QapQpf3^ +C {Qpqf Qa/3 «, /3 = 1 . . . 3. 

(35) 

We use standard results from the theory of elliptic partial differential equa- 
tions to deduce that Q* is actually a classical solution of the system ([35]) and 
Q* is smooth and analytic everywhere in [TO]. Given smooth boundary 
conditions, Q* is also smooth up to the boundary. 

3.2.2 Upper bounds for the order parameters 

The global minimizer Q* can be expressed in terms of a pair of eigenvectors 
(n*,m*) and the scalar order parameters (s*, r*), as in ([5]), 

Q* = s* (n* n* - -I ) + r* ( m* m*- -I ) (36) 



3 / V 3 
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and 



\Q*\^ = l(^s*^ + r*^ -s*r*) . (37) 



We partition the (s, r)-plane into three regions: (a) Ri = {{s, r)\s,r > 0} - 
the top quadrant, (b) R2 = {{s, r)\s < 0; r > s} and (c) R3 = {r <0; r < s}. 
In we have the inequahties 

-(s*+r*)2 < |Q*|2 < -(s*+r*)2. (38) 
6 3 

Similarly, for R2, we have that 

-(r* -2s*f <\Q*\^ <-(r* -2s*f (39) 
6 3 

and for R3, 

- (s* - 2r*f < |Q*|2 < - (s* - 2r*f . (40) 
6 3 

For every 77 > 0, we define the bounded region = {(s*, r*) ; |Q*| < 
in the (s, r)-plane. Let be the isosceles triangle in the (s, r)-plane with its 
vertices at the points - {(r/, 0) , (0, ry) , (— r/, —r])}. Then it follows immediately 
from daHl), (I39D and (00]) that 

^.fl.^^v^^VE, (41) 



so that C necessarily implies that C T^. Our first result in this 
section is an explicit upper bound for the norm of a global energy minimizer 
in the low-temperature regime a < ^ and as the preceding discussion 
shows, this upper bound allows us to define the admissible domain for the 
equilibrium scalar order parameters. 

Theorem 1. Let Q* ^ A be a global minimizer for the energy functional 
Ilg, where A and Ilg have been defined in (T^ and [W\) respectively. We 
work in the temperature regime a < ^ and make the following assumption 
about the boundary condition Qq, 



I 6 + - 24ac 1 

IQ»WI<"-| .Til (^2) 

where |Q| has been defined in [T^). The condition is equivalent to 
requiring that the boundary order parameters are contained inside and 
the bulk triangles A(r) defined in Section 3.1. Then Q* obeys the following 
global upper bound on Q - 



mr)\<'-±^^^^ /orrea (43) 
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Proof. We assume that the contrary holds i.e. the subset Q* ■ 

has positive measure. The subset fi* clearly does not intersect dil. since 
the boundary condition Qo obeys this upper bound from assumption 
We define a perturbation, Q, of the global minimizer Q* as follows 



re 17; |Q*(r)| > ^^+.^1^} c 



Q(r) 



'Q*(r), r€n\n*, 



where 



We note from (ETI) that 



b + Vb"^- 24ac 
2V6c 



(44) 



(45) 



(46) 



by definition. It is evident from (1440 that Q agrees with Q* everywhere 
outside Qs* and hence belongs to our admissible space. Moreover, Q has 
constant norm on the set fi* i.e. |Q(r)| = F for r G 17*. 
We obtain an upper bound for the free energy difference 

/lg[Q]-/lg[Q*] = / /b(Q) + L|VQ|2-(/b(Q*) + L|VQ*|2) dV (47) 

JO.* 



where fs is as in ([T 

We can explicitly compute VQ 

r 



as shown below 



VQ(r 



Q*(r) 



VQ 



*|2 



IQ 



7^ (Q;,Qkfc) {^kk^kk) < |VQ*(r 



(48) 



smce 



< 1 on 17* by definition. 



|Q*(r) 

Consider the function G : [0,oo) 



defined by 



G{u) 



-u 



a 

2 ~ 3^ 



u + -u' 
4 



(49) 



We estimate the bulk energy density difference in terms of the function G 
as follows 



fsiQ) - fB (Q*) = -trQ' - -trQ3 + - trQ 



2 ^ 



IQ 



*|2^ 



2 3 
6trQ*3 
3TQ~ 



|trQ*^ - hrQ*' + | (trQ*^; 



* 13 



r3 _ |Q*|3) + J (r^ - |Q*|^) < G(|Q*|) -G{T). 



(50) 
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In the last step of ([50]) . we use the equahty 

trQ^ = r^TT^, 

|Q*|3 ' 

(r^ - |Q*P) < on W and the inequaUty 

W - 71' ^''^ 

from Lemma [H One can readily verify that G{u) attains a local maximum 
for u = T and G'{u) < for all u>T. Therefore, 

G(|Q*|)-G(r) <0, (52) 

since |Q*| > F on fi* by definition. 

We substitute (|l8]), ([50]) and ([52]) into ^ to obtain 



Ilg 



Q -/lg[Q*] <0, (53) 



contradicting the absolute energy minimality of Q* . We thus conclude that 
ri* is empty and 

|q-(r)|<r=^±^^^ (54) 

for all points r G 0. □ 
Lemma 1. Let /3(Q) he defined as follows - 

(irQ2 

Then < /3(Q) < 1. 



/3(Q) = 1-6^7^ Qg5o. (55) 



Proof. The quantity /3(Q) is known as the biaxiality parameter in the liquid 
crystal literature and it is well-known that /3(Q) G [0, 1] [15]. We present a 
simple proof here for completeness. 

Since 6 ^ 3 > 0, the inequality /3(Q) < 1 is trivial. To show /3(Q) > 0, 

we use the representation ([8]) to express trQ'^ and trQ^ in terms of the order 
parameters s and r. 

trQ^ = ^ (2s^ + 2r^ - 3sV - Ssr^) 
9 

trQ2 = ^ + - sr) (56) 

A straightforward calculation shows that 

(trQ^)^ = — (45^ + - 12s V - 12sr^ + 26s - SsV^ - 3s V^) 
81 
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and 

(trQ2)^ = A (56 ^ ^6 _ 3<,5^ _ 3^^5 _ 7<,3^3 ^ gg2^4 ^ gg4^2^ _ 

One can then directly verify that 

(trQ2)^ - 6 (trQ3)^ = 2s {s - rf > (57) 

as required. □ 

As a further illustration, let us assume that there exists an uniaxial global 
energy minimizer in the admissible space, A in (|14p . where the boundary 
condition Qo is of the form 

Qo = so (^no no - ^I^ , no : 50 ^ 5^ (58) 

and < So < min{s-|_, 1} is a positive constant. Then we have 

Lemma 2. Let Q„ be an uniaxial global minimizer of Ilg ii^- the admissible 
space A, with a smooth, uniaxial and physically realistic boundary condition 
Qo, as in 158\) . Then Q„ is necessarily of the form 

Qu = Su ^n„ (g) n„ - -il^ (59) 

for some function s„ : 17 ^ M and unit-vector field n„ : 17 — > 5"^. The 
equilibrium scalar order parameter is non-negative everywhere and obeys the 
following inequalities 

< Su{r) < s+ r G n. (60) 

Proof. We prove Lemma [2] by contradiction. Let 17* = {r G 17; Su(r) < 0} 
be a measurable interior subset of 17. The subset 17* does not intersect 917 
by virtue of our choice of Qo in ([55]) . Consider the perturbation 



Quir) 



Q„(r), rGl7\17*, 

(61) 

^-Q„, r G 17*. 



Then Q^ G A and Qu coincides with Qu everywhere outside 17*. We ex- 
plicitly estimate the free energy difference, Ilg[Qu] — IlgIQu], as shown 
below 

IlgIQu] - IlgIQu] = [ fsiQu) + L\VQu\^ - [fe {Qu) + ^IVQJ^) dV = 

Jn* 

= [^trQldV= [^sldV<0 (62) 
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since trQ^ = |s^, < on 17* by assumption and 6 > 0. This contradicts 
the absolute energy minimality of Q^. Hence, il.* is empty and > 
everywhere in 0. 

The upper bound in (j60p follows directly from (j43p i.e. 



\Qu\ = Y ^ Y ^•^+ 

and since s„ > 0, we have that < < s+ from (|63p . Lemma [2] now 
follows. □ 

Remark: We are not guaranteed the existence of a uniaxial global en- 
ergy minimizer in our admissible space. We assume the existence of in 
Lemma [2j A more technically precise formulation of the problem would be 
to minimize Ilg in the restricted class of uniaxial Q-tensors 

X = |q G {n- So); Q = s (^n n - ^I^ a.e. in 

where s is a real-valued function and n is a unit-vector field, subject to 
uniaxial boundary conditions. In this case, one can prove the existence of a 
uniaxial global minimizer Qu in the restricted class Au and the statement 
of Lemma El still holds. However, the proof is technically more involved and 
we omit the details for brevity. 

Given the explicit upper bound F in (j43p for the norm of a global en- 
ergy minimizer, the corresponding equilibrium scalar order parameters are 
confined to the bounded region Sr = {(s, r); |Q*| < L} in the {s, r)- 
plane. We define the elastic triangle Aez(T) to be the triangle Tygp, where 
T C Sr C in (|4ip . The elastic triangle can be explicitly specified in 

terms of the temperature and the material-dependent bulk constants. One 
can directly verify that Sr C Aez(T) C if and only if 

6 - c 6^ 

- — -+T* <T <^— + T*. (64) 
6a 24ac 

and Sr D T D if and only if 

V 2^ 

T <^(h-2c)+T*. (65) 
3a 

In other words, for temperatures T G [Oi 3^ ~ 2c) +T*), the equilib- 
rium order parameters may move outside the physical triangle. For the 
liquid crystal material MBBA, Sr D for T < 44.52°C i.e. the Landau- 
De Gennes predictions fail to be consistent with the probabilistic second- 
moment definition of Q within a 2°C-neighbourhood of the nematic-isotropic 
transition temperature. 
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3.3 Maximum principle approach 

In this section, we carry out a parallel analysis in the high-temperature 
regime a > ^ and extend our analysis to more general Landau-De Gennes 
energy functionals, where fs is a general even polynomial in the Q-tensor 
components. 

Theorem 2. Let Q* be a global minimizer of Ilg the admissible class 
A, where Ilg o,nd A have been defined in U^) and (T^ respectively, in the 
temperature regime « > Then the function |Q*| : 17 ^ M attains its 

maximum on the domain boundary. In particular, if the boundary condition 
Qo satisfies 

|Qo(r)|<-^ re 90, (66) 

then the scalar order parameters of Q* are contained inside T^, for this 
high-temperature regime. 

Proof. The proof proceeds by contradiction. We work in the temperature 
regime a > ^ and assume that the function |Q*| : O ^ M attains a strict 
maximum at an interior point r* G Q, where |Q*(r*)| > 0. 

The global minimizer Q* is a classical smooth solution of the Euler- 
Lagrange equations 

2LAQa/3 = aQa/3+b (^t^Q'^Sa/S - QapQpfS^ +C {Qpqf Qa/3 L > 0, a, (3=1... 3. 

(67) 

Therefore, the function |Q*P : O ^ M is also a smooth function and we 
necessarily have that 

A|Q*|2 < at r* G 0, (68) 

according to our hypothesis 0. We compute A|Q*p at this interior maxi- 
mum point. One can readily show that 

LA|Q*|2 = 2L (|VQ*|2 + Q^Q*^. . (69) 



We substitute the Euler-Lagrange equations (|67j) into (j69j) to obtain the 
following 

LA|Q*|2 = 2L|VQ*|2 + (a|Q*|2 - 6tr (Q*)^ + c|Q*|^) + ^ (trQ*^) Q*^ 

= 2L|VQ*p + atrQ*^ - 6trQ*^ + c (trQ*^)^ , (70) 

since Q*j = 0. 

Consider the function G : 5o — > M defined by 

2\2 



G(Q) =atrQ2-6trQ3 + c(trQ2)\ (71) 
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Then G is bounded from below by 

G(Q) > /i(IQI) = a|Q|2 - ^IQI' + clQl" 

from Lemma [TJ The function /i : — > M has its global minimum at the 
isotropic state, Q = 0, and 

M|Q|)>0 Q/0 

in the temperature regime o > This implies that 

G(Q*) > at r* G 

and consequently A|Q*p(r*) > from ()70p . This contradicts our hypothesis 
and Theorem [2] now follows. 

In particular, if the boundary condition Qo satisfies the hypothesis (|66p . 
then the global energy minimizer Q* satisfies the inequality 

IQ*(r)l < ^ (72) 

on 0. From (j41|) . this is sufficient to ensure that Q* is physically realistic, in 
the sense that its scalar order parameters do not take values outside T^. □ 

Theorem [5] shows that in the high temperature regime a > ^ , the norm 
of a global energy minimizer attains its maximum on the boundary. The 
isotropic state Q = is the global minimizer of the bulk energy density, 
fs, in this high-temperature regime and it is not surprising that we observe 
a dissipation of order in the interior. However, it is interesting that there 
are no local fluctuations in the interior i.e. there are no interior regions 
where |Q*| experiences a local increase compared to the boundary norm - 
maxrgon |Qo(r)|. Therefore, Theorem [2] suggests a monotonic decrease in 
order as we move away from the boundary and it would be interesting to 
analytically estimate the characteristic length scale of order decay for this 
problem. 

Our methods readily extend to a more general bulk energy density, fB,n, 
which is a polynomial of even degree 'n' in the Q-tensor invariants i.e. trQ^ 
and trQ^ [16j with n > 4 (since fB,n has to be minimally quartic to allow a 
first-order nematic- isotropic phase transition). We take fB,n to be 

fB,n = a2(r)trQ2-a3trQ3+a4 (trQ2)^+. . .+ ^ a„,p (trQ^)™ (trQ^) 

m,p£Z+;2m+3p=n 

(73) 

where Z^" denotes the set of non- negative integers, 03, 04 > and 

af,o > Yl \am,p\- (74) 

m,p£Z+ ;p> 1 ;2m+3p=n 
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The first coefficient 02 (T) lias a linear dependence on the absolute tempera- 
ture by analogy with (fT9|) whereas the remaining coefficients {03, 04, . . . , {am,p}} 
are taken to be temperature-independent, material-dependent bulk con- 
stants. 

We define the corresponding Landau-De Gennes energy functional to be 

In[Q] = [ /b,„(Q) +L|VQp dV (75) 
Jn 

and our admissible space is 

An = {Qe W^''' {Vt- 5*0) ; Q = Qo on dVt] (76) 
where the Sobolev space is defined to be [7] 

VF^'"(n;S'o) = |q G 5o; y IQI" + |VQr < ooj (77) 

and Qo is a smooth, physically realistic boundary condition in the sense of 
()66p . We have the following result by analogy with Theorem [H 

Proposition 2. Let Q* he a global minimizer of In in the admissible space 
An- Then 

|Q*| < max < C (a2, as, 04, . . . , {a^ „}) , maxIQol/' on fl (78) 
where C is a positive constant that only depends on the absolute temperature 



and the bulk coefficients and is independent of the elastic constant L in { 75). 



Proof. The existence of a global energy minimizer, Q*, for /„ in the space 
An, follows from the direct methods in the calculus of variations and we 
can use standard results in elliptic regularity to deduce that Q* is smooth 
everywhere in Q. Consider the function |Q*P : 0, — > R'^ and assume that it 
attains its maximum at the interior point r* G 0,. 

The global minimizer Q* is a classical solution of the Euler-Lagrange 
equations 

orrk dfB,n ^dfB,nj- 

'^^^^•'^'^ - aq- - SdQ^J'^ ^^'^ 

where the second term is a Lagrange multiplier accounting for tracelessness. 
We multiply both sides of 1^9^ by Qij and use A|Qp = 2 (| VQp + QijQij,kk) 
to get 

LA|Q|2 = Q,,.|^ + 2L|VQ|2. (80) 



Then 

\*|2 



A|Q*r < at r* G (81) 
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from our hypothesis. 

We note that Qtj is a polynomial of degree 'n' in Q and from 
Lemma [Hand (I74p . we have that 



> A-(|Q| 



where 



i^(IQI) = «2|Q|'-^|Q|' + ...+ an^o- ^ |am,p| I IQI 

'm,p£Z+;p>l;2m+3p=n 



(82) 

The function ii' : ^ R is a polynomial of degree 'n', has 'n' zeros 
{|Qi|) IQnl) • • • ) IQn|} where |Qi| < IQ2I < • • • < |Qn| and K is a monoton- 
ically increasing function of |Q| for |Q| > |Q„|. 

If |Q*(r*)| > |Q„|, then we necessarily have that A|Q*p > at r* (from 
([80|) ). contradicting the hypothesis ([HT]) . We, thus, conclude that 



IQ1 < IQnl on n (83) 

where |Q„| = C (02, as, 04, ... , {am,p}) can be explicitly expressed in terms 
of the bulk coefficients. Proposition [2] now follows from combining (j83p and 
the maximum norm of Qo on the boundary. □ 

The explicit upper bound (1780 allows us to define the admissible domain 
for the equilibrium scalar order parameters as in Section 13.21 For certain 
choices of the bulk coefficients, this domain is larger than and conse- 
quently, the equilibrium scalar order parameters may take values outside 
the physical triangle. 

4 Discussion 

We have studied qualitative properties of global minimizers of the Landau- 
De Gennes energy functional, /lGi in smooth three-dimensional geometries 
with Dirichlet boundary conditions. We have obtained an explicit upper 
bound for the norm of a global energy minimizer in terms of the temperature 
and material-dependent bulk constants, independent of the elastic constant. 
In particular, we have defined two triangles in the order-parameter, (s, r)- 
plane: (a) the bulk triangle A(T) which accounts for the stationary points 
of /b (b) the elastic triangle Aei{T) which accounts for the effects of the 
elastic energy density and A(T) C Ae/(T). The equilibrium scalar order 
parameters take values inside or on the boundary of A(,i(T) and the distance 
D{T) between /\f,i{T) and A(r) scales as 



D{T) <s+ = . (84) 
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This, in effect, quantifies the effect of elastic perturbations on the bulk 
energy minimizers. Secondly, this explicit bound is also compared to the 
probabilistic bounds in ([7|) and we find that the equilibrium scalar order 
parameters may move outside the physical triangle and take physically 
unrealistic values larger than unity, in the low-temperature regime. For the 
liquid crystal material MBBA, the Landau-De Gennes predictions fail to 
be physically realistic within a 2°C-neighbourhood of the nematic-isotropic 
transition temperature. 

A natural question is - how can we reconcile the differences between the 
Landau-De Gennes predictions and the probabilistic second-moment defi- 
nition of Q in the low-temperature regime? The Landau-De Gennes bulk 
energy density, fs, has no term that enforces the probabilistic bounds in ([7]) 
or penalizes configurations that lie outside T^. A first-step in this direction 
is to use a Ginzburg-Landau approach [3] . We define a modified Landau-De 
Gennes energy functional, as shown below 

Fe[Q]= [ fB{Cl) + fe{m+L\VCl\^dV (85) 



where 



0, IQI<^, 



(86) 



Jb is as in (fT7|l and e > is a small positive parameter. We can obtain 
an explicit upper bound for the norm of a global energy minimizer using a 
maximum principle approach as shown below. 

Proposition 3. Let Q* be a global minimizer of in the admissible space 
A in |T^ . Then 



1 be^ JU + 16e2 (c - 6a) + (62 _ 24ac) 

Q < max < —= , —= 1 ■= , max 

\^/% \/6(8 + 2e2c) ^/6(8 + 2e2c) rean ' 

(87) 

where a,b,c are the bulk constants in fs in 

Proof. The existence of a global energy minimizer Q* follows from the direct 
methods in the calculus of variations and we use standard arguments in 
elliptic regularity to deduce that Q* is smooth everywhere in 0, up to the 
boundary. Then |Q*P : 0, — > M"*" is a smooth function and we assume that 
it attains its maximum at an interior point r* G 17 and |Q*(r*)| > (i.e. 
/, / at r* from ([86])). It follows that 

A|Q*p < air* £ n. 
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The global minimizer Q* is a classical solution of the Euler-Lagrange 
equations 



Repeating the same arguments as in the proof of Theorem [2l we have the 
following inequality 

LAIQp > M(IQI) = a|Qp - -^|Q|=^ + c|Q|4 + ^ f IQ'^ ^ 



We study the function M : O — > M above. This function has precisely 
three zeros: 

IQI = 



_ fee^ y^64 + 16e2 (c - 6a) + (62 _ 24ac) 

' ' ~ ^/6(8 + 2e2c) ^ ^/6(8 + 2e2c) ^^^^ 

a„d M(IQI) > fo. IQI > ^ + Therefore, 

we must have 



\Q*(r*)\ < , V64 + 16e2 (c - 6a) + (62 _ 24^ 

' ^' " ^(8 + 2e2c) ^(8 + 2e2c) 

in order to have A|Q*|2 < at r* G 17. 
We, thus, conclude that 



|Q.| < 1 , ^ "'^ + V64+16.Me-6,)+.Mi,^-24a.) 

[\/6 ^/6(8 + 2e2c) ^/6(8 + 2e2c) redn' 

(90) 

where the second term accounts for the maximum of |Qo| on d^l. □ 

In the limit e ^ and for a physically realistic boundary condition Qq 
satisfying the upper bound ([57|) reduces to 

|Q*|<^ + 0(e) on 11 (91) 
V6 

and consequently, the equilibrium scalar order parameters take values within 
for all temperature regimes. Further, one can show that the stationary 
points of the modified bulk energy density 

/B,e(Q) = /B(Q)+/e(|Q|) 

are either isotropic or uniaxial Q-tensors as in Proposition [T] and fB,e also 
predicts a first-order nematic-isotropic phase transition. Therefore, the mod- 
ified Landau-De Gennes energy functional reproduces all the qualitative 
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features of Ilg in (fT6]) . whilst respecting the probabihstic bounds ([7]) in the 
Umit e ^ 0^, for all temperature regimes. 

However, the Ginzburg-Landau approach in ()85p does not contain any 
information about the probabilistic second- moment definition of Q in ([3]). 
A more systematic approach is given in [11] where they define a modified 
bulk energy density from the Maier-Saupe free energy Ims in (|11P - 

^b(Q)= inf / ^/'(p)logV'(p) dp-ic/(r)|Q|2 (92) 

where 

= |V' G (5'; M+) ; Q = p - ll^ V(P) c^p} (93) 

is the space of all probability distribution functions ip that have a fixed 
normalized second moment Q, as in (l3|0- 

The first term in (j92p is the entropy contribution, where we minimize the 
integral over all probability distribution functions that have a fixed second 
moment Q. This term diverges whenever the probabilistic bounds in d?]) are 
violated and enforces the equilibrium scalar order parameters to lie strictly 
inside T^. The second term in (|92p is simply the Maier-Saupe interaction 
energy. The uniaxial case is treated in [TT] . We plan to study the biaxial case 
and include spatial inhomogeneities into this model. This will be reported 
in future work [2j. 
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